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PLOT SUBROUTINE 
FOR FORTRAN 
WITH FORMAT 

(1620-FO-004 Version 1) 



Author: Edward A. VanSchaik 
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DECK KEY 

1. FORTRAN Plot Demonstration Deck to Read 
Values from and Compute the Sine, Cosine 
and Natural Log 

2. Plot Subroutine Program Deck 

3. Plot Subroutine Source Deck 



4. Plot Subroutine Demonstration Data Deck 



Plot Subroutine for FORTRAN with Format 
(1620-FO-004 Version 1) 



Author: Edward A. VanSchaik 
IBM Corporation 
370 West First Street 
Da5rton 2, Ohio 

Direct Inquiries To: Author 



A. Purpose /Description: This program presents a subroutine which can be 
easily used in FORTRAN for plotting single or multiple curves without 
the GOTRAN problem of determining which curve has the lower magnitude . 
Also this method can be used to generate the curves on cards for off-line 
printing. 

B. Method : N/A 

C. Restrictions and Range : Normally a range of 0-70 should be the limitation 
for typewriter output since the carriage can handle only 80 to 85 positions. 
The range for card output should be 72 characters minus the digits in the 
independent variable. 

D. Accuracy : Negative numbers, floated numbers, and numbers outside of 
the allowable range produce an ERROR F9 message. . 

E. Machine Conficruration ; Basic 1620 Card system with minimum of 20K 
memory. The program is designed to work with the FORTRAN system 
without automatic divide. 

F. Program Requirements: The subroutine requires 432 core locations. 

G. Source Language : Fortran 

H. Program Execution Time: N/A 

I. Check Out Status : N/A 

J. Sample Problem Running Time : Approximately 45 minutes. 

K. Comments : Major modifications must be made to use this subroutine 
with the FORTRAN version with automatic divide. This program and 
Its documentation were written by an IBM employee. It was developed 
for a specific purpose and submitted for general distribution to interested 
parties In the hope that It might prove helpful to other members of the 
. data processing community. The program and Its documentation are 
essentially in the author's original form. IBM serves as the distribution 
agency in supplying this program. Questions concerning the use of the 
program should be directed to the author's attention. 




Plot Subroutine of IBM 1620 FORTRAN with Format 



INTRODUCTION 

Since many IBM 1620 applications are in areas where a rough plot of the 
results can be quite useful, it was felt that a subroutine incorporated into 
FORTRAN was necessary. GOTRAN provides a statement for plotting 
requiring a separate program and a knowledge of GOTRAN language. This 
write-up presents a subroutine which can be easily used in FORTRAN for 
plotting single or multiple curves without the GOTRAN problem of deter- 
mining which curve has the lower magnitude. Also this method can be used 
to generate the curves on cards for off-line printing. 

METHOD OF USE 

Since a method of plotting multiple curves with any special characters was 
desired, the normal method of output (PRINT or PUNCH) was utilized with 
the incorporation of a new type FORMAT character Z. This subroutine 
is used with the output routines rather than the normal arithmetic statement. 
The basic FORTRAN statements to use this plot subroutine are: 

PRINT 2, X, M, N, J 

2 FORMAT (FIO. 5. 3Z ) 

The first variable in the list is the independent variable which can be either 
fixed or floated point number. This variable will be handled in the regular 
output routine where the value is printed. The Z character in the format 
accomplishes the following; 

1. The remaining values in the list will be plotted instead of printed. 

2. The Z is handled by the compiler in the same routines designed for 
the H character except that the address compiled will be the address 
in the symbol table corresponding to the subroutine whereby the 
FORTRAN input/output package will branch to this address then to 
the relocated subroutine. The number and type of plot characters 
are stored similar to the H type. 

3. The remaining values in the list are assumed to be fixed variables 
previously scaled to the range of output record minus the digits in 

the independent variable. (Example: for the case shown, the numbers 
should be scaled between and 87 minus 10 (or 77). Normally a 
range of 0-70 should be the limitation for typewriter output since the 
carriage can handle only 80 to 85 positions. The range for card 
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output should be 72 characters minus the digits in the independent 
variable. Negative numbers, floated numbers, and numbers 
outside of the allowable range produce an ERROR F9 message. 

4. The first character after the Z will be used to plot the first 
variable in the list to be plotted. Remaining list values will use 
the next character until the list is satisfied. Therefore, no 
limitation on the number of variables exists except the 72 char- 
acter statement size. If the number of variables in the list 
exceeds the number of plot symbols, an ERROR F9 message will 
occur. 

5. The range of the plot starts at the next position after the last digit 
of the independent variable. Therefore, if the value of the variable 
to be plotted is 0, the plot symbol will be inserted in the first 
position after the independent variable. 

The statement may be either PRINT or PUNCH depending on whether the 
output is on the typewriter or cards. The subroutine sets up the FORTRAN 
output area with the independent variable and all plot symbols before giving 
the write command. For typewriter output, a record mark is inserted in 
the next position after the highest plot value to conserve typing time 
(A minimum of 20 characters is always typed). This produces a much 
faster method than GOTRAN and eliminates all problems of relative 
magnitude of curves. 

FORTRAN Processor Modifications 

To make the plot routine useable on a 20K card system, space was used 
where the ACCEPT TAPE and PUNCH TAPE statements were previously 
located. If these commands are called, an error message indicating a 
missing format number will occur. These changes are: 

Card No. 3024 

Location New Instruction Card Col. 



2030 
2042 
2054 



CM 

BNE 

B 



14 17323 00069 
47 04022 01200 
49 02174 



21-32 
33-44 
45-52 



Location 



New Instmction 



Card Col. 



2162 
2174 
2186 
2198 
2210 



B 

TDM 
TPM 
TR 
B 



49 02218 00000 

15 08719 00001 

16 03493 19870 
31 17322 17324 
49 04118 



13-24 
25-36 
37-48 
49-60 
61-67 



Card No. 3053 
Location 



New Instruction 



Card Col. 



4070 
4082 
4094 
4106 



TFM 
BNE 
TR 
TFM 



16 04243 00031 
47 02030 01200 
31 17322 17324 
16 0349 



31-42 
43-54 
55-66 
67-72 



Card No. 3054 



Location 



New Instruction 



Card Col. 



4112 
4118 



TFM 



3 04264 
16 03620 04174 



3-8 
9-20 



These patches allow the processor to check for the Z format, insert the 
subroutine name location, and set for subroutine load. The instructions 
at 2174 and 2186 assume that the subroutine will be the first added after 
the normal six from IBM. If this is not true, then the addresses of 2174 
(8719 for subroutine load) and 2186 (T9870 for subroutine address) must 
be changed. The normal rules for adding a subroutine must be followed. 

Plot Subroutine 

Enclosed at the end of this write-up is an SPS listing and flow charts of 
the Plot Subroutine for the FORTRAN subroutine deck without automatic 
divide. For use with automatic divide, the address for the symbols of 
the input /output package must be changed. This includes CON2 which refers 
to the instruction in SWC. 

The subroutine requires 432 core locations. The trailer card requires 
a value of 00434 in Col.l- 5 as stated in the FORTRAN write-up. 

General Comments 



With format control, the output plot Ccin be quite complete. The following 
technique which can be used in FORTRAN provides a means of printing 
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the graph scales: 



PRINT 5, 0, 50, 100, 150, etc. 
PRINT 5, 1, 1, 1, 1, 1, 1, etc. 
5 FORMAT (I 11, 15, 15, 15, etc.) 

FORTRAN will print the variables in the list as the value of each quantity. 
The 111 can be made so that the first value prints after the independent 
variable. Example (using previous example) 

. 50 100 150 200 250 300 350 400 
111111111 
. 50000 * . + 

Although the previous sections of this write-up indicates that only one 
independent variable may be printed, this is not true. Assume it is 
desired to print x and y and also plot y. This could be done as follows: 

PRINT 6, X, Y, MY 
6 FORMAT (FIO. 5, FIO. 5, IZ*) 

This will print as follows: 

X y 
+XXX. xxxxx+xxx. xxxxx * 

Where X and Y are printed and MY (Y after scaling) is plotted. The range 
of MY must be reduced to keep the over record with limits . The general 
rule is that any variable to the left of the first plot variable in the list will 
be printed under normal FORTRAN rules. 

The name PLOT must not be used in the program since the compiler will 
assume this to be a function name. 

Modifications for Extra Memory 

To convert this program for machines with 40K or 60K, the compiler 
modifications previously described must be altered. The digit in column 
44 of card 3026 (shown as 1 ) must be changed to 3 or 5 depending on the 
amount of additional memory. 
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Sample Program 



Included with this write-up is a sample program and listing to illustrate 
the use of this subroutine. This program works as follows using data 
values of 1 to 89: 

1. Switch 1 on - switch 2 off 

SIN, COS and LGN are printed on the typewriter and 
punched on cards. i 

2. Switch 1 off - stops typewriter output. 

3. Switch 1 on - switch 2 on 

SIN, COS and LGN are pimched on cards. 
SIN and COS are plotted. 

It should be noted that without typewriter output, this program permits 
almost complete overlap of input, output, and computing providing an 
effective demonstration of the 1620 card buffering. 
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S83 00 C 
583 00 C 
58300 C 
58300 C 
58300 100 
583 20 1 
58340 11 
58364 
58508 
58652 
58660 12 
58684 20 
58708 
58720 
58744 
58780 
58828 
58852 
58876 
58900 



S AMPLE PROBLEM LISTING 

1620 F0RTRAN DEM0NSTRAT I0N T0 READ VALUES FR0M 1 T0 89 

AND C0MPUTE THE SINE,C0SINE AND NATURAL L0G 0F THESE NUMBERS 

SWITCH 1 0N T0 0BTAIN TYPEWRITER 0UTPUT 

SWITCH 2 0FF T0 LIST VALUES 0N T0 PL0T SIN AND C0S 

IF(SENSE SWITCH 1) 10,20 

IF(SENSE SWITCH 2) 11,12 

PRINT 6 

PRINT 5,0,1,2,3,4,5,6,7,8,9,10 
PRINT 5,1,1,1,1,1,1,1,1,1,1,1 

G0 T0 20 
PRINT 1 
PUNCH 1 

00 150 N=1,89,2 

READ 2,1 

A=l 

X=A*3. 1415926 /I80. 
B=SINF(X) 
C=C0SF(X) 
D=L0GF(A) 
PUNCH 3,1 ,B,C,D 
IF(SENSE SWITCH 1) 110,150 
IF(SENSE SWITCH 2) 120,130 
MPL0T = 50.* B 
NPL0T = 50. *C 
PRINT 7,1 ,MPL0T,NPL0T 
G0 T0 150 
PRINT 3,1 ,B,C,D 
C0NTINUE 
PRINT 4 • 
PAUSE 
G0 T0 100 

F0RMAT{3H I 6X6HSIN(I) 6X6HC0S(I) 6X6HLGN(I)) 
F0RMAT(I5) 

F0RMAT(I5,F12.7,F12.7,F12.7) 

F0RMAT(37HRUN C0MPLETE , PRESS START F0R NEXT SET ) 

F0RMAT(I7,I5,I5,I5,I5,I5.I5,I5,I5,I5,I5,I5,I5 ) 

F0RMAT(24HBASE=I ,*=SIN( l),+=C0S( I)) 
F0RMAT( I6,2Z*+) 
END 



58980 110 
59000 120 
59048 
59096 
59144 
59152 130 
59212 150 
59248 
59272 
59284 
59292 1 
59436 
59458 
59496 
59594 
59676 
59748 
59782 

PR06 SW 1 0NF0R SYMB0L TABLE, PUSH START 

J9999 SIN 

J9989 SINF 

29979 C0S 

39969 C0SF 

59959 ATAN 

[9949 ATANF 

[9939 EXP 

J9929 EXPF 

r9919 L0G 

J9909 L0GF 

J9899 SQRT 

[9889 SQRTF 

[9879 PL0T 

59869 PL0TF 

J9859 5100 

J 9849 5010 

J9839 5020 

J9829 5011 

I198I9 5012 .6 
f9809 5006 
J9799 5006 
T9789 5005 



J9779 5005 

19769 S005 SAMPLE PROBLEM LISTING (Con't . ) 

3759 DOOl 
,9749 0002 
J 973 9 5005 
[9729 5004 
J 97 19 5003 
I9709 5oo5 
29699 5007 
39689 5008 
29679 5002 
39669 5010 
39659 5001 
39649 5001 
3963 9 5 1 50 
39629 N 
29619 S002 
39609 0002 
39599 I 
29589 A 
39579 X 

39569 5141592651 
29559 SOO 
39549 T800000053 
39539 B 
29529 C 
39519 D 
59509 5003 
59499 5003 
39489 5110 
39479 5120 
59469 5130 
39459 MPL0T 
59449 '5000000052 
59439 NPL0T 
19429 5007 
[9419 5007 
59409 5004 
J9399 5004 

SW 1 0FF T0 IGN0RE SUBR0UTINES, PUSH START 
PR0CESSING C0MPLETE 



SAMPLE PROBLEM OUTPUT 



L0AD DATA 



I 
1 
3 
5 
7 
9 
11 
13 
15 
17 
19 
21 
23 
25 
27 
29 
31 
33 
35 
37 
39 
^1 
k3 
kS 
^7 
k9 
51 
53 
55 
57 
59 
61 
63 
65 
67 
69 
71 
73 
75 
77 

Z9 
81 

83 

85 

§7 
89 



SIN{I) 
.017'^524 
.0523359 
.0871557 
.1218693 
.156434if 
. 1 908089 
.22if9510 
.2588190 
.2923716 
.3255681 
.3583679 
.3907311 
.4226182 
.4539904 
.4848096 
.5150380 
.5446390 
.5735764 
.6018150 
.6293203 
.6560590 
.6819983 
.7071067 
.73 13 536 
.7547095 
.7771459 
.7986355 
.8191520 
.8386705 
.8571672 
.8746196 
.8910064 
.9063077 
.9205048 
.9335804 
.9455185 
.9563047 
.9659258 
.97^3700 
.9816271 
.9876883 
.9925461 
.9961946 
.9986295 
9998476 



C0S(l) 
.9998476 
.9986295 
.9961946 
.9925461 
.9876883 
.9816271 
.9743700 
.9659258 
.9563047 
.9455185 
.9335804 
.9205048 
.9063077 
.8910065 
.8746197 
.8571673 
.8386705 
.8191520 
.7986355 
.7771459 
.7547095 
.7313537 
.7071067 
.6819983 
.6560590 
.6293204 
.6018150 
.5735764 
.5446390 
.5150381 
.4848096 
.4539906 
.4226182 
.3907312 
.3 583680 
.3255682 
.2923717 
.2588190 
.2249510 
.1908090 
.1564345 
.1218694 
.0871558 
.0523360 
.0174524 



LGNd) 

.0000000 
1.0986122 
1.6094379 
1.9459101 
2.1972245 
2.3978952 
2.5649493 
2.7080502 
2.8332133 
2.9444389 
3.0445224 
3.1354942 
3. 21887 58 
3.2958368 
3.3672958 
3.4339872 
3.4965075 
3.5553480 
3.6109179 
3.6635616 
3.7135720 
3.7612001 
3.8066624 
3.8501476 
3.8918202 
3.9318256 
3.9702919 
4.0073331 
4.0430512 
4.0775374 
4.1108738 
4.1431347 
4.17't3872 
4.2046926 
4.2341065 
4.2626798 
4.2904594 
4.3174881 
4.3438054 
4.3694478 
4.3944491 
4.4188406 
4.4426512 
4.4659081 
4.4886363 



RUN C0MPLETE, PRESS START F0R NEXT SET 
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SAMPLE PROBLEM OUTPUT 



BASE=l,*=SIN(l),+=C0S(l) 

1 2 3 4 5 6 
1111111 
1* 

3 * 

5 * 

7 * 

9 * 
11 * 
13 * 
15 * 
17 * 
19 * 
21 * 
23 * 
25 * 
27 * 
29 * 
31 * 
33 * 
35 * 
37 * 
39 * 
41 
43 
45 
47 
49 

51 + 

53 + 

55 + 

57 . + 

59 + 

61 + 

63 + 

65 + 

67 + 

69 + 

71 + 

73 + 

75 + 

77 + 

79 + 

81 + 

83 + 

85 + 

87 + 

89+ 

RUN C0MPLETE, PRESS START F0R NEXT SET 



* + 

+ 
+ * 




Store No. of 
Plot Symbols 
in T3+2 



Plot Subroutine 
Flow Diagrams 

First Entry into 
Subroutine by 
Z Character. 



Locat e First j 
Plot Symbol 
(Plot-f 35) 



Alter Swc. 
Instr uction to 
Allow Re-entr|r 
to Plot 1 



Reduce Count 
of Remaining 
Output Area (F ) 
By 2 



Set'Ov^put to i 
Insert at i 
Min. 6f QZ-irA) 




Return to Main 
Program through 
FORTRAN Output 
. Routine . 
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Re- Entry to Subroutine 
for each plot symbol. 




FROM PREVIOUS 




PAGE 



Store Plot 
1 Symbol in 

1 Output Area 




Restore 
SWC 

Instruction 

. 



Branch to Write 
Output Section of 
FORTRAN output. 
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ERROR MESSAGE 
ROUTINE 
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SOURCE LISTING 



00 RG 




05000 


PL0T TF 


*+35,BLK10+11 ,0 


05000 26 05035 07275 






AM 


*+23,2,OlO 


05012 Tl 05035 OOOIS2 






TF 


T3+2 


05024 26 07428 00000 
05036 Tl 05035 000U2 






AM 


PL0T+35,2,O1O 


TF 


SWC+6,C0N1+6,1 


05048 2S 04998 05420 


AM 


BLK5f6,2,10 


O5O6O 11 06326 OOOU2 






SM 


P,2,10 








05072 12 07807 000ri2 


TFM 


PT+l8,QZ+40 




0^84 16 07778 n80<)^ 


B 


BRBK 


050Q6 49 07360 00000 


D0RG 




05104 


PL0T1 SM 


T3+2,1 ,10 


0<*104 12 07428 OOOS1 




BN 




05116 1f7 05378 01300 




TR 


RX-9,Z0TZ-1O 


05128 3 t' 00089 07994 




TF 


*+23,SWC-l ,0 


051 'to 56 05163 04991 




TF 


RX 


05152 26 00098 00000 


BNF 


ERR0R,RX-3,O 






05164 ^4 05378 00095 


c 


RX.P 








05176 24 00098 07807 
05188 ^6 05378 01300 


BNL 


ERR0R,,O 


BNF 


*+24,RX,0 








05200 ^4 05224 00098 
05212 Ti9 05378 00000 


B 


ERR0R,,O 


TF 


PL0T2+6,BLK5+6,O 


05224 26 05278 06326 




A 


PL0T2+6,RX,O 




05236 2l 05278 00098 


A 


PL0T2+6,RX,O 








05248 1\ 05278 00098 


TF 


PL0T2+11,PL0T+35,O1 




05260 IE 05283 05035 
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SOURCE LISTING 



PL0T2 


TF 


0,0 












05272 26 00000 


00000 




C 


PT+18,PL0T2+6,1 








*+24,,0 


05284 25 07778 






BNL 






PT+18,PL0T2+6,1 


05296 <C6 053 2n 






TF 










05308 25 07778 


05278 




60 


PL0T3,T2,O 






05320 53 05364 


05365 




TFM 


SWH+6,BRBK 








05332 16 07850 


57360 




TF 


SWC+6,C0N2,1 










05344 2S 04998 


0543 2 




B 


SWD 




D0RG 




05356 49 04692 
05364 


00000 




*-3 




PL0T3 


AM 


PL0T+35,2,O1O 










053 A4 Tl ncn^c 


UUUUZ 




BB 














053 7 A 4? nnnnn 
05378 






D0RG 


*-9 




ERR0R 


RCTY 
WATY 


ERC0NT 


05378 34 OOOOn 

/ w \J\J\J \J\J 


00102 






PL0T2+48,,O 


05390 39 07977 


00100 




B 




ERC0NT 






05402 ^9 0^20 


00000 


OS 


,7977 










07977 00000 




T2 


OS 


,5365 










05365 00000 




Z0TZ 


OS 


,8004 




C0N1 




08004 00000 




B 


PL0T1 ,,0 






C0N2 






05414 ?9 05104 


00000 


DC 


7,^307876 








,7264 


05432 00007 




BLKIO 


DS 










07264 00000 




T3 


OS 


,7^26 










07426 00000 




SWC 


OS 


,^992 










04992 00000 
00098 00000 




RX 


OS 


.98 




P 


DS 


,7807 


07807 00000 
06320 00000 












BLK5 


DS 


.63 20 
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SOURCE LISTING 



PT 


DS ■ 


.7760 


QZ 


DS 


,8053 


BRBK 


DS 


,7360 


SWD 


DS 


,4692 


SWH 


DS 






DEND 





END 0F PASS I I 



07760 00000 
08053 00000 
07360 00000 
04692 00000 
076kk 00000 
00000 
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